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Abstract 

We consider the effect of compactification of extra dimensions on the onset of classical chaotic "Mix- 
master" behavior during cosmic contraction. Assuming a universe that is well-approximated as a 
four-dimensional Friedmann-Robertson- Walker model (with negligible Kaluza-Klein excitations) 
when the contraction phase begins, we identify compactifications that allow a smooth contraction 
and delay the onset of chaos until arbitrarily close to the big crunch. These compactifications are 
defined by the de Rham cohomology (Betti numbers) and Killing vectors of the compactification 
manifold. We find compactifications that control chaos in vacuum Einstein gravity, as well as in 
string theories with ]\[ = 1 supersymmetry and M-theory. In models where chaos is controlled 
in this way, the universe can remain homogeneous and flat until it enters the quantum gravity 
regime. At this point, the classical equations leading to chaotic behavior can no longer be trusted, 
and quantum effects may allow a smooth approach to the big crunch and transition into a sub- 
sequent expanding phase. Our results may be useful for constructing cosmological models with 
contracting phases, such as the ekpyrotic/cyclic and pre-big bang models. 
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I. INTRODUCTION 



The behavior of spacetime near a big crunch singularity has been a topic of research for 
many decades. The classic studies of Belinskii, Khalatnikov, and Lifshitz (BKL) □ D □ 
and others have shown that the contraction to the crunch can 

proceed either smoothly or chaotically. Chaos arises when the universe is unstable to small 
inhomogeneities and anisotropies in curvature or matter fields. These "dangerous" pertur- 
bations grow and eventually dominate the dynamics, driving the universe to an anisotropic 
state, expanding along some axes and contracting along others. The axes and their rates 
of contraction jump to new values when new curvature or matter terms grow to dominate. 
The jumps generally repeat an infinite number of times before the big crunch itself. The 
chaotic, oscillatory evolution to the big crunch is often known as "Mixmaster" behavior or 
"BKL oscillations." 

Models that generalize four dimensional Einstein gravity have been classified by whether 
they necessarily exhibit chaotic behavior near a big crunch. This classification is established 
assuming "generic" initial conditions, in which there is a finite, but possibly small, energy 
density in all fields present in the model. Additionally, one assumes that the classical Einstein 
equations remain valid up to the big crunch itself. Under these assumptions, it has been 
shown that vacuum Einstein gravity in spacetime dimension less than eleven, as well as all 
uncompactified ten dimensional string theories and M-theory, will inevitably suffer chaos as 
the big crunch is approached jj, 0, ^| . So too will Einstein gravity containing only perfect 
fluids with equation of state w < 1, where w = p/p is the ratio of the fluid's pressure to its 
energy density. On the other hand, there are cases where chaos is not inevitable. Among 
these are Einstein gravity with a free massless scalar field, leading to w — 1, and a universe 
containing a perfect fluid with equation of state w > 1 [1^ . 

The presence of chaos during gravitational collapse is a potential problem for cosmological 
models with a big crunch /big bang transition, such as the ekpyrotic/cyclic and pre big- 
bang scenarios 3, Q, till Q • In these models, the universe undergoes collapse to a big 
crunch, followed by a transition to the conventional big bang and subsequent expansion. 
It is assumed that, during the collapsing phase, the universe is nearly homogeneous and 
isotropic, with a scale invariant perturbation spectrum. BKL oscillations arising during the 
collapsing phase would destroy homogeneity and isotropy, producing a chaotic spacetime 
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with structure down to arbitrarily small scales. In this situation, a big crunch/big bang 
transition is unlikely to be describable in a deterministic manner, and it is questionable 
whether a homogeneous and isotropic universe with the long range correlations required by 
observations could emerge jl9j . Thus, avoiding chaos is an essential feature of cosmological 
models with a collapsing phase. 

We take the point of view that avoiding chaos all the way to the the big crunch is too 
restrictive a condition for viable cosmological models. One expects classical general relativity 
to break down at a small but finite time before the big crunch is reached, perhaps of order a 
Planck time tpL or string time ts- After this point, quantum effects become significant and 
we can no longer trust the classical physics that predicts a chaotic approach to the big crunch. 
Provided the universe evolves smoothly and non-chaotically until tpi, it is conceivable that 
quantum gravity effects allow the universe to pass smoothly through the big crunch and into 
a subsequent expanding phase. For example, recent work [20( has revealed that the degrees 
of freedom present in string and M-theory (extendend objects such as strings and branes) 
can evolve smoothly through certain types of big crunch singularities. This suggests that 
chaotic behavior is absent in the quantum gravity regime, and furthermore a nonsingular 
transition from a big crunch to a big bang is possible. 

Our focus in the present work is on the classical evolution of the universe before tpL, and 
whether it is possible for it to evolve smoothly so long as we can trust the classical equations 
of motion. The evolution of the universe through the subsequent quantum regime, and the 
transition to an expanding phase, are important though unsettled issues. In this work we 
have nothing to add on these topics. However, if chaos can be controlled during the classical 
evolution of the universe, there remains the hope that the subsequent quantum evolution will 
preserve the long range correlations, isotropy and homogeneity so essential for cosmology. 

In this paper, we consider models that are well described by a classical, four-dimensional 
effective field theory long before the big crunch. The four-dimensional metric is that of 
a nearly homogeneous and isotropic Friedmann-Robertson-Walker (FRW) universe, with 
small perturbations to the metric, matter and Kaluza-Klein fields, and with the Hubble 
radius H" 1 much larger than the compactification length scale R c . We study the evolution 
of chaotic behavior as the universe contracts, including the effects of all massive Kaluza- 
Klein modes, and thus all of the degrees of freedom of the higher-dimensional theory. We 
show that the emergence of chaos can be controlled provided the compactification manifold 
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M satisfies certain topological conditions. For these topologies, the dangerous perturbations 
that formerly led to chaos acquire masses of order the inverse of the compactification length 
scale, m ~ 1/R C - The presence of mass terms slows the growth of energy density in these 
fields, and prevents them from becoming cosmologically relevant so long as the Hubble 
parameter is larger than their mass. When the time until the big crunch becomes less than 
R c , or equivalently m < H, the suppression ceases to operate, and the energy density in 
the dangerous modes can grow at their usual unsuppressed rate. However, since the energy 
density in these heavy modes has been greatly suppressed relative to light modes up to this 
point, they cannot dominate the energy density until the universe has contracted further 
by an exponential factor. Typically, the massive modes do not dominate before we enter 
the quantum gravity regime at roughly tpL, at which point the classical evolution equations 
cannot be trusted. In these circumstances, we say that chaos has been "controlled." 

In this paper, we focus on a classical effect which reduces the importance of chaos in 
compactified models. This is especially relevant for models, such as string- or M-theory, in 
which the compactification of extra dimensions is an essential element. An excellent example 
is given by eleven-dimensional supergravity, whose bosonic sector contains a four-form field 
strength in addition to the metric. Without the four-form, pure eleven-dimensional gravity 
is not chaotic. For some choices of the topology of the compactification manifold, it is 
possible to remove the light modes of the four-form field. For these topologies, the previously 
chaotic eleven-dimensional supergravity theory will behave like the non-chaotic eleven- 
dimensional pure gravity theory. 

More generally, when one studies a fully uncompactified model, one finds that chaos arises 
from dangerous modes that are nearly spatially homogeneous along all dimensions. The en- 
ergy density in these modes scales rapidly enough to dominate the universe and cause chaos. 
As we detail below, for some choices of the compactification manifold, the classical equations 
of motion forbid spatially homogeneous modes along the compact directions, for topological 
reasons. In the four-dimensional effective theory, this is reflected in the appearance of large 
mass terms for the associated degrees of freedom. As we will show, as long as the four- 
dimensional Hubble parameter is larger than their mass, the energy density in these massive 
modes grows far more slowly than the energy density in the light modes which dominate the 
dynamics. As the Hubble radius falls below the compactification length scale, the energy 
in the massive modes begins to grow more quickly, but due to its relative suppression, it 
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remains dynamically irrelevant all the way to the Planck time. In the present work, we focus 
exclusively on the classical evolution of fields and find that it suffices to control chaos. It 
is possible that quantization, by imposing a further constraint on the initial perturbations, 
would further suppress chaos. The quantum production of heavy Kaluza-Klein modes is, 
naively at least, completely negligible all the way to the Planck time. 

In this work, we consider both pure Einstein gravity and models with additional matter 
fields. We focus on p-form matter fields, with exponential couplings to a scalar "dilaton" 
field 0, defined by the action, 

s = ~ T^Ti)! / e ^ Fp2+1 ^ d ° x ' Fp+1 = dApi (1) 

with A a constant [21]. Supergravity and string models commonly include p-form fields 
with couplings of this type. In the following, we will always use "p" to denote the number 
of indices on the gauge potential A p . While many models that generalize four dimensional 
Einstein gravity contain fermionic fields, throughout this work we will focus exclusively on 
the bosonic sector. We will also neglect more exotic terms in the p-form action such as 
Chapline-Manton couplings and Chern-Simons terms, which at any rate we do not expect 
to be relevant for chaos Q. Another important type of matter, the perfect nnid, has been 
discussed elsewhere [14], and is not affected by the compactification of extra dimensions. 

In Section |Hj we review some results regarding the emergence of chaos during gravita- 
tional collapse that are required in later sections. This section is primarily concerned with 
distinguishing between chaotic and non-chaotic models. Most importantly, we define the 
gravitational, electric and magnetic stability conditions that must be satisfied if a theory 
is to avoid chaotic behavior. In Section IIH[ we introduce some key aspects of models with 
"controlled chaos," that are the subject of the current work. We establish that giving masses 
to dangerous modes prevents them from causing chaos. We then describe how the masses of 
these fields are determined by the compactification manifold. With these results established, 
we present our central new results in Section HVl These are the "selection rules" that deter- 
mine a subset of stability conditions that must be satisfied in order to control chaos. These 
rules express the precise correspondence between the properties of the compactification man- 
ifold, and the chaotic behavior of the lower-dimensional theory after compactification. In 
Section we give examples with specific compactification manifolds that are able to control 
chaos in vacuum Einstein gravity, string theories with M = 1 supersymmetry, and M-theory. 
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We summarize our conclusions in Section fVTl and suggest some areas for further research. 



II. REVIEW OF GRAVITATIONAL AND p-FORM CHAOS 

The essential principle underlying the emergence of chaos is that, near a big crunch, 
solutions to the Einstein equations are strongly unstable to perturbations. This phenomenon 
may be studied using a suitably general metric, such as the generalized Kasner metric, 

D-1 

ds 2 = -dt 2 + J2t 2pj (°j) 2 , (2) 

3=1 

where D is the dimension of spacetime, the <jj = aj fl (x)dx fl are independent of time, and the 
big crunch occurs at t = 0. The Kasner exponents pj may be spatially varying, but upon 
substituting (J2J) into the Einstein equations, one finds the pj are constrained by the Kasner 
conditions, 

D-1 D-1 

£Pi = 1 > 5>? = L ( 3 ) 

3=1 3=1 

The first condition defines the Kasner plane, the second defines the Kasner sphere, and we 
may term their intersection the Kasner circle. For these anisotropic metrics it will be con- 
venient to define the analogue of the Hubble parameter in the conventional, isotropic FRW 
solution. Using (J2J), the metric on equal time hypersurfaces is given by h^ v = . sh- 
allowing us to define, 



H — — log Vh, where h = det h^. (4) 

(JjL 

The first Kasner condition implies that H = t^ 1 for any choice of the Kasner exponents. 
We will find this Hubble parameter a useful guide to the typical dynamical timescale of the 
gravitational field. 

The Kasner metric (J2J) has been widely used as a tool to understand the behavior of 
"generic" spacetimes near a big crunch singularity , 0, Q . It is an 

approximation, to leading order in t, of an exact solution of the Einstein equations, in 
which we have neglected the influence of spatial derivatives and curvature terms. To check 
whether this approximation is consistent, we substitute the Kasner metric into the Einstein 
equations, and check that terms corresponding to spatial derivatives and curvature appear at 
subleading order in t. This corresponds to spatial derivatives and curvature terms becoming 



irrelevant as the big crunch is approached. It has been shown rigorously 

7] 

and numerically 
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12j that, when these terms are subleading, solutions to Einstein's equations asymptotically 
approach Kasner form as t — > 0. 

A useful feature of generalized Kasner universes is that the conditions determining 
whether the curvature terms are irrelevant can be expressed entirely in terms of the Kasner 
exponents. The Einstein tensor for the Kasner metric (0) may be split into purely time 
derivative terms, and terms arising from the curvature of the spatial slices. Details of this 
decomposition are given in One finds that the time derivative terms all scale as 

t~ 2 , while the terms arising from spatial curvature scale as t~ 2<KPi+P3 ~ Pk \ for all triplets k. 
Dangerous components of the curvature are those whose corresponding terms in the Einstein 
equations grow more rapidly than t~ 2 , as t — > 0, thus invalidating the Kasner approxima- 
tion. These dangerous curvature terms are absent provided that the gravitational stability 
conditions, 

Pi + Pj —pk<l, all triples k, (5) 

are satisfied. If the stability conditions are satisfied, then the evolution is guaranteed to be 
smooth and Kasner-like all the way to the big crunch. These conditions turn out to be very 
restrictive; in the absence of matter, it is only possible to simultaneously satisfy (JHJ) and (jnj) 
when the spacetime dimension is greater than ten j^J. 

When the gravitational stability conditions are not satisfied, then spacetime will exhibit 
chaotic behavior. Violation of the gravitational stability conditions (0) means that the 
generalized Kasner solution (J2J) is invalid, and so we must find a different description. One 
useful picture recasts the evolution of the metric in terms of geodesic motion of a point mass, 
undergoing reflections from a set of sharp walls 8]. The free flight of the point between wall 
collisions is described by the Kasner metric, where the Kasner exponents give the momentum 
components of the moving point. Collisions with walls correspond to spatial curvature terms 
temporarily dominating the Einstein equations, and result in a sudden change in the Kasner 
exponents. The motion of the point in the chamber defined by these walls is chaotic, and 
so the dynamics of spacetime is as well. 

It is useful to distinguish more precisely between the various possibilities with respect 
to the stability conditions A model is chaotic if the stability conditions cannot be 
satisfied for any choice of Kasner exponents satisfying the Kasner conditions. We will also 
term models chaotic when the stability conditions are only satisfied at isolated points on the 



Kasner circle. For example, it is possible to almost satisfy the stability conditions in any 
spacetime dimension with the so called "Milne" solutions, in which a single Kasner exponent 
is unity and the rest are zero. This leads to marginally dangerous curvature components 
that scale exactly as t~ 2 . Thus, one might argue that if these curvature components are not 
dominant initially, they will remain subdominant all the way to the crunch, and chaos will 
not arise. However, since this occurs only for isolated points on the Kasner circle, any small 
perturbation of the Kasner exponents away from the Milne solution results in a violation of 
the stability conditions and the emergence of chaos. These solutions are thus not practically 
useful from the perspective of avoiding chaos in cosmological models, since they do not admit 
the small inhomogeneities and anisotropics that must be present in any physically realistic 
scenario. Therefore, we consider a model non chaotic only when there exists an open region 
of the Kasner circle in which all of the stability conditions are satisfied. 

The presence of matter can either enhance or suppress chaos. Three important examples 
are a free massless scalar field, p-form fields, and perfect fluids. The first, a scalar field, 
suppresses chaos by modifying the Kasner conditions. A homogeneous, massless scalar field 
coupled to the Kasner metric will evolve as, 

<P(f) = 0o +£Vlogt, (6) 

with the constants 0o and p^ determined by the initial conditions. Including the stress 
energy from in the Einstein equations results in the new Kasner conditions, 

j=l 3=1 

While the scalar field "Kasner exponent" p^ enters the Kasner conditions, it does not enter 
the gravitational stability conditions It is now possible to find pj that satisfy these 
stability conditions in any spacetime dimension D. For example, the isotropic choice, 

ID -2 1 

satisfies all of the stability conditions. Moreover, there is a finite, open neighborhood on the 
Kasner circle surrounding the isotropic solution where the stability conditions are satisfied. 
Essentially, there are two key properties of the scalar field that enable it to suppress chaos. 
The scalar field has an isotropic stress energy tensor, and therefore does not enhance any 
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preexisting anisotropy in the Kasner metric. Also, the scalar field energy density scales as 
t~ 2 , and thus grows sufficiently rapidly to remain relevant near the big crunch. 

The addition of p-form fields coupled to tends to enhance chaos. A p-form field strength 
with p > 1 has an anisotropic stress energy tensor, which tends to enhance any preexisting 
anisotropy in the Kasner metric. If the p-form field dominates the energy density, then 
it tends to drive the universe to anisotropic oscillations and then chaos. A homogeneous 
p-form evolves simply, and its dynamics depend on whether it has a time index (an electric 
p-form) or whether all indices are spatial (a magnetic p-form). Using the equation of motion 
and the Bianchi identity, one finds, 

(F P+ J hh ~ jp - 9r ~ const - (9) 

Using these solutions, we may construct the stress energy tensor for the p- form field, and 
compare its time dependence with the t~ 2 leading dependence of the homogeneous terms 
in the Einstein equations. Note that, unlike the scalar field, we do not include the gravita- 
tional back-reaction from the p-form fields, and merely check if it remains subdominant if 
subdominant initially. One finds that the p-form energy scales more slowly than t~ 2 , and 
therefore cannot cause chaos, when the following p-form stability conditions are satisfied; 

^2Pj-^Y>° (electric), (10a) 

p 

y^^Pj — — — < 1 (magnetic). (10b) 
p+i 

The electric conditions involve a sum of p distinct Kasner exponents, corresponding to the p 
spatial indices of an electric p-form field strength. The magnetic conditions likewise involve 
a sum over p + 1 distinct Kasner exponents. Thus, for each p- or (p + l)-tuple of Kasner 
exponents, there will be a corresponding stability condition. The inequalities (fTUj) are often 
referred to individually as the electric or magnetic stability conditions. If they are violated, 
then chaos will arise. 

The inclusion of a matter component with w > 1 can suppress chaos ^| . The w > 1 fluid 
suppresses chaos by growing rapidly to dominate the energy density of the universe. Once it 
dominates, the power-law Kasner solution (0) is replaced by a solution of the approximate 
form, 

3=1 
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where the Cj are arbitrary constants. When w > 1, this solution converges rapidly to isotropy 
as t — > 0. Near the big crunch, this solution may be thought of as a Kasner universe with 
Kasner exponents, 

2 

Pi = (D-l)(l + wY (12) 
The gravitational stability conditions are clearly satisfied in this case, and thus chaos does 
not arise. The p-form case is somewhat more subtle, and depends on whether we realize the 
w > 1 fluid by introducing a potential for the dilaton 0, or as a separate matter component 
that does not couple to the p-forms. In either case, for any p-form and coupling A, there is 
a corresponding w cr i t (p, A), for which chaos is eliminated when w > w crit . 



III. MASSIVE MODES, CHAOS, AND COMPACTIFICATION 

In this section, we will distinguish a subclass of the chaotic models, those with "con- 
trolled" chaos. The salient feature of controlled chaos is that dangerous modes are sup- 
pressed, relative to t~ 2 , for an adjustable epoch of cosmic history. In these models, dangerous 
modes acquire a mass m, and their contribution to the total energy density of the universe 
is suppressed so long as m/H > 1. These masses arise through compactification, and are of 
order -R" 1 , with R c the characteristic length scale of the compactification manifold M.. In 
a universe with Kasner metric 0, H = £ _1 . Thus, in these models chaos cannot emerge 
until t <C R c , and is is possible to ensure that chaos does not arise before the Planck time 
tpL- In Section UlI Al we discuss how the growth of dangerous modes is suppressed by mass 
terms, and estimate the suppression factor. We give a simple argument based on a scalar 
field in an isotropic universe, leaving a discussion of the general p-form case to Appendix 
1. Sections IIII Bl and IIII CI describe the mechanism by which compactification gives the re- 
quired masses to p-form and metric degrees of freedom. In both cases, massless modes in the 
lower dimension can only exist when M. admits p-form or vector fields with certain special 
properties. In Section IIV1 we use these results to give the "selection rules" that express 
the correspondence between these special properties of M. and the chaotic properties of the 
compactified theory. 

In the following, we assume that spacetime has the form R > 3 x A4, with Ai a compact 
manifold. Indices M, N, P . . . denote directions in the total spacetime R 1 ' 3 x Ai, while 
/i, v... denotes directions along R 1,3 and m,n,p... along Ai. Gmn denotes the total 
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metric, with h^ v the metric on R 1,3 and f mn the metric on jVl. In later sections we will 
need to distinguish between coordinate indices and tangent space indices on R 1,3 x M.. We 
use A, B, C... for tangent space indices along the full space, a, j3... for those along R 1,3 , and 
a,b... for those along M.. 



A. Massive and massless modes 



The suppression of massive modes in a collapsing universe may be illustrated using the 
equation of state w. This is defined as the ratio of the pressure to energy density for a 
perfect fluid, 



P T j 



w = - = ^ ( 13 ) 
where we assume that we are in the comoving frame where the stress energy tensor is 
diagonal. As we are primarily interested in cosmological models, it is sufficient to consider the 
case where the universe is an isotropic, four dimensional FRW model after compactification. 
Conservation of stress energy implies that the energy density p of a perfect fluid with equation 
of state w depends on the scale factor a as, 

p = p a-^ 1+w \ (14) 

where po is the energy density when a is unity. In a contracting universe, the component 
with the largest w grows most rapidly, and eventually dominates the total energy density. 
A homogeneous, massive scalar field has a perfect fluid stress energy with equation of state, 

4> 2 -m 2 4> 2 

w = - . (15) 

<p 2 + m 2 (j) 2 

From equations ((H)) and ()15j) it is readily seen that the energy density of a massive scalar 
field must scale more slowly than that of a massless one. A massless scalar field will always 
have w — 1, and thus its energy density p scales with a as p ~ a" 6 . The energy density of a 
massive scalar will scale with an effective w between zero and unity, depending on the ratio 
m/H. When m/H ^> 1, far from the big crunch, the scalar field's dynamics is dominated 
by the mass term in its potential. Using (•) to denote the time average, the virial theorem 
implies that (<fi 2 ) = m 2 ((j) 2 ), and therefore w = [2^. Thus the energy density in the 
massive field scales as p ~ a -3 , far more slowly than the massless field. Near the crunch, 
when m/H 1, the mass term has a negligible effect on the field's dynamics. In this limit, 
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tpL tj 

logt 



to 



FIG. 1: Scaling of the energy density in dangerous modes before and after compactification. With- 
out dangerous mode scales more rapidly than t 2 and eventually dominates the energy 
density of the universe, as illustrated in A. With a mass, the dangerous mode scales much more 
slowly than t~ 2 , until a time t\ at which m = H, as seen in B. After this point, the mode will 
scale as usual. The chaos is "controlled" when the dangerous energy density cannot catch up to 
1/t 2 before the Planck regime is reached at tpL- 

<fi 2 ^> m 2 <f) 2 , and w approaches unity from below. In this regime, the energy density in 
massive and massless fields will scale identically with time. 

Using this equation of state argument, we can estimate the exponential suppression of 
the energy density in massive fields. The basic process by which dangerous modes are 
suppressed, and then grow near the big crunch is illustrated in Figure Q Although we have 
discussed only the scalar case above, the same w = scaling when m/H > 1 obtains for 
general p-form fields, as discussed in Appendix 1, and so this estimate applies to those 
fields as well. It is important to emphasize that, while the energy density in massive modes 
is always growing, it grows more slowly than the t~ 2 scaling required for the field to be 
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cosmologically relevant. Therefore, if we wish to estimate the importance of a given energy 
component, we should consider the ratio of the energy density in the component to t~ 2 . This 
quantity, t 2 p(t), may be thought of as measuring the ratio of the energy density in a given 
component to the total density, for, 

where we have used Planck units. Only when t 2 p(t) is increasing as t —>■ can a given 
component grow to dominate the universe. 

We begin by choosing a reference time t , a t the beginning of the contracting phase. 
We consider a model where the four-dimensional effective theory begins to contract with a 
background equation of state w, so that, 

a(t) = (t/t ) W, (17) 

where we have normalized a = 1 when t = t . Combining this with (|14|). one finds that for 
massive modes, 

t 2 p(t) = [t 2 p(t)] Q - (t/t )^, (m/H > 1), (18) 

where [t 2 p(t)} denotes t 2 p(t) evaluated at t . This equation is valid up until m/H ~ 1, 
when the mass terms are becoming irrelevant. Denoting by t\ the time at which m = H, 
one finds, 

T = RcH , (19) 
where we have taken m = 1/R C . We then have, 

[t 2 p(t)] 1 = [t 2 p(t)] Q ■ (R c Ho)^, (m/H ~ 1). (20) 

This equation shows the suppression in the fractional energy density in massive modes 
during the period in which m/H > 1. The suppression is controlled by the ratio of the 
compactification length scale R c to the Hubble horizon Lh = l/-ffo at the beginning of the 
contracting phase. 

As we approach the big crunch, we have m/H < 1, and the dangerous modes can grow 
as usual. We parameterize this growth by an exponent 6, so that, 

t 2 p(t) = (t/t 1 )- s [t 2 p(t)] 1 , (m/H<l). (21) 
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when 5 > 0, a mode will grow and eventually dominate the energy density of the universe. 
One can see by comparing (}2"T]) to our discussion in Section |H] that 5 is merely twice the 
amount by which a given mode violates the stability conditions. Typically, S will be of order 
one. Now we define a time t eq , at which the dangerous modes have grown sufficiently so 
that the fractional energy density in dangerous modes is equal to that at the beginning of 
the contracting phase, or, 

[t 2 p{t)] eq =[e P {t%. (22) 

Using ipn jl and (f2Tjl . we find, 

t eq = t {R c H ) 1+ W+^ (23) 

Finally, we define a time t^ om , at which the dangerous modes formally dominate the universe, 
corresponding to t 2 p(t) ~ 1. Then one finds, 

tcU > m = t eq ([t 2 p(t)} ) 1/5 . (24) 

Chaos is controlled provided that the dangerous modes do not dominate before a Planck 
time from the big crunch, corresponding to td om <tpL- 

Having established the formulae we will need for our estimate, we may now insert reason- 
able values for our variables. Let us assume the contraction phase begins when the Hubble 
parameter is of order the present value Hq (as occurs in ekpyrotic and cyclic models). Then, 
1/Hq ~ 10 61 Lpl. As a first example, we assume that w = 1 during the contraction as 
is characteristic of compactifications of Kasner universes. A typical value for S that arises 
from working with the p-form spectrum of string models and gravity is 5 = 2, although the 
precise value of 5 will of course depend on the specific model under consideration. If we take 
R c ~ 10 • Lpi, as an example, then we find the suppression factor (j2U|) at H ~ m to be, 

(R C H )^ ~ 1(T 60 . (25) 

The dangerous modes grow to have the same fractional energy density as they had at the 
beginning of the contracting phase at the time, 

t eq ~ 10~ 29 t PL . (26) 

Thus, the dangerous modes cannot grow to be even as relevant as at the beginning of 
cosmic contraction, until the universe is well within the quantum regime. We would need 
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to approach even closer to the big crunch for these modes to dominate the universe, but at 
this point we no longer expect our classical equations to be valid. 

As another example, we consider a case with large compact dimensions, where m is of 
order the weak scale, 10~ 16 Mp^, which corresponds to taking R c ~ 10 16 Lp^. Now we find 
that the suppression when the Hubble radius equals the compactification scale is, 

(R C H )^ ~ lO" 45 , (27) 

and, 

t eg = 10" 6 t PL . (28) 

Again, we will be well within the quantum gravity regime before the dangerous modes can 
potentially dominate. If «) > 1 during the contraction phase, as occurs in ekpyrotic and 
cyclic models, the dangerous modes are suppressed by a much greater factor than in these 
examples, as is evident from the expressions above. 

As a final relevant issue, note that we have taken the mass m of the dangerous modes 
to be constant in time. Generally we expect that the mass m of a given field will be time 
dependent. This occurs since the field's mass is determined by the compactification manifold 
A4, and M. will be evolving with time during cosmic contraction. Let us say that the mass 
m(t) evolves as, 



m{t) = m ^-j , (29) 

during the contracting phase, where for simplicity we take w = 1. Then, because the energy 
density will go like p ~ m/a 3 , we find that 

*V=[*U(rY* ( 30 ) 



However, the time t\ at which m/H ~ 1 is now, 



t{ +b = X (31) 
m 

and therefore the suppression when m/H ~ 1 is, 

j^ji = (H /mo) ~ (H R C , ), (32) 

where R c ^ is the characteristic length of Ai when t = t . This is precisely the factor found 
in the case where m is constant in time ([20)1 . The difference is that the time ti, at which 
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m/H ~ 1, shifts. This shift compensates for the different growth rate of p, with the net 
result that the suppression factor remains the same. 

From the equations above, it is clear that there is a potential problem if b < — 1, which 
would lead to t 2 p increasing during the contracting phase. However, our assumption that 
the universe is of Kasner type excludes this possibility. Field masses m are related to the 
compactification length scale R c by m ~ 1/R C - In a Kasner universe, we expect that R c ~ t p , 
with p a Kasner exponent, or average of several Kasner exponents. However, the Kasner 
conditions (f3*|) or (|7j) imply that p < 1, and thus m must vary more slowly than 1/t. This 
implies that b > —1, and so the suppression operates as before. 



B. The p form spectrum 



Having established that the massless modes are the only relevant modes near the big 
crunch, we now describe how these massless modes are determined in terms of the compact- 
ification manifold A4. The story is familiar from the study of higher-dimensional models of 



particle physics 
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24. 



251 ] although we discuss some special features of the time dependent 



situation which must be taken into account. An useful feature of the p-form mass spectrum 
is that the existence of massless modes is determined entirely by the topology of A4, and 
not by its metric. This simplifies the task of finding manifolds M. that lead to controlled 
chaos, since we need only specify their topological properties. 



1. Time independent compactification 

First, we will review the situation for the time independent case. For clarity, we will 
neglect here the exponential coupling to the dilaton field, which amounts to an overall 
multiplication by e A * of the Lagrangian density. The coupling is fully accounted for in the 
analysis given in Appendix 1. The action for a p-form gauge potential A p in 4+n dimensions 
is, 

S = - 1 „ [ dA p ■ dA p V^G d 4+n x } (33) 
{p+l)lj 

where A p can depend on all coordinates, and has indices along both R 1,3 and Ai, 

A v = \AJx^,x m )] . (34) 
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The conventional compactification analysis begins with an expansion of the p-form A p as, 

A ?= E £^ )A & (i) (35) 



r+s=p i 



using a basis of s-forms /9s , with indices and coordinate dependence only along Ai, 

P? = W(x m )] aia2 ... as . (36) 

The abstract index i labels the s-form under consideration, and when Ai is compact it 
takes discrete values, infinite in number. The "coefficients" in this expansion are r-forms 
a?', that depend on the noncompact coordinates on R 1,3 , and have indices along R 1 ' 3 only, 

a f = . (37) 

It is convenient to choose the gauge d^ = 0, and select the fijfi to be eigenfunctions of 
the Hodge-de Rham Laplacian A = dd^ + d^d on Ai, with eigenvalues \i , 

A/?« = A»fl«>. (38) 

The are normalized so that, 

/ (3i l) ■ (3^ V^f d n x = 6 ir (39) 
Jm 

Substituting the expansion (J33)) into the original action (JBlIj) results in, 

S = ~J2 1 1V [ " ^ + X s )a r ] " "9 ^ ( 40 ) 

i,s ( r ^ ''J 

(i) 

where we have rescaled the a r hy a constant in order to canonically normalize the kinetic 
terms. This demonstrates that a single p-form in 4 + n dimensions yields many r-forms ar 
after compactification, whose masses mft are related to the eigenvalues A^ by [rnf 1 ^ = \$ . 
These are the "Kaluza-Klein" modes. The operator A has a positive semi-definite spectrum 
on manifolds which, like Ai, have a Euclidean metric. Therefore the effective masses are all 
real. 

The (p — s)-forms with zero effective mass are determined entirely by the topology of 
Ai, and not by its metric structure. As discussed above, massless (p — s)-forms arise from 
s-forms (3 S satisfying = 0, conventionally termed "harmonic" forms. Hodge's theorem 



26, 



271 ] states that the number of harmonic s-forms is equal to the dimension of 7i s {At), the 
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s th de Rham cohomology class of Ai. The dimension dim7Y s (Ai) is also known as the s th 
Betti number of Ai, conventionally denoted b s (Ai). This is a topological invariant, which 



does not change under smoot 



i deformations of Ai and its associated metric structure. The 



Poincare duality theorem 



28, 



291 ] gives a simple geometric interpretation of these cohomology 



classes; the quantity dim?-^ s (Ai) counts the number of s-dimensional submanifolds that can 
"wrap" Ai, and cannot be smoothly contracted to zero. In this counting, two submanifolds 
are considered equivalent if one can be smoothly deformed into the other. Thus, for every 
inequivalent noncontractible submanifold of Ai with dimension s, a p-form gives rise to a 
massless (p — s)— form field after compactification. 



2. Time dependent compactification 

The case where the compactification manifold Ai changes with time introduces new 
features, but in the end does not substantially modify the conclusions reached above. The 
main difference is that it is no longer possible to assume that the eigenbasis of forms (3^ 
defined by (|35j) depends only on the compact coordinates x m . In particular, the forms will 
depend on time. This introduces additional cross terms which must be taken into account. 
Below, we will neglect the variation of Ai along directions x other than time. We denote 
by d\_M the exterior derivative tangent to the manifold Ai. Thus, 

= d\ M Pi <) + dtAP®. (41) 

We may now use our freedom to choose the basis modes j3s \ and define the modes Ps(t) with 
X s l \t) > to be the instantaneous eigenforms of the Hodge-de Rham operator, restricted 
to act on Ai only, 

A\ M p®(t) = \®(t)$\t). (42) 

We find it convenient to relax the requirement that the zero modes Ps(t) with As = be 
eigenforms of the Hodge-de Rham Laplacian. Instead, we will merely require that they be 
representatives of the de Rham cohomology of Ai. Inspection of the reduced action shows 
that this condition is sufficient to guarantee that the zero modes still result in massless form 
fields. Furthermore, we adopt the normalization convention, 

/ PP-PP V^fd n x = 5 ij f y^fd n x. (43) 

J M J 
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This differs from the usual normalization convention ()39|) by only a multiplicative constant 
in the static case. It has the advantage of not introducing any spurious time dependence of 
the /3a from the changing volume of M.. Maintaining this normalization condition requires, 

f $®-(3i j) V^fd n x = 0. (44) 
Jm 

With these conventions and definitions, we find that the p-form action (}3*3*j) splits into two 
parts, Si and S 2 , with 

Sl = ~ 2 / da r ] ■ da r ] + X s\t)a? ■ d 4+n x (45) 

s,i ^ ''J 

and, 

S 2 = -J2 T^W / A #°) • («(«% A #°) + • ■ #°) ^ rf 4+ "x, (46) 

s,i 

where terms that are identically zero due to mismatched indices are not included. Again, 
we have rescaled the ai to obtain canonically normalized kinetic terms. The terms in S2 
threaten to substantially modify the action in the time dependent case. However, these 
terms vanish or are negligible. The first term is zero due to our normalization convention 
(J43|) and its consequence (|44|). The second term is a contribution to the effective mass of 
ai l \ For the A« = modes, the representatives /3g of the de Rham cohomology are time- 
independent, and so the /3 2 terms vanish. When As > 0, the additional contribution to the 
effective mass will be positive, but since these modes are already massive it will not change 
the qualitative features of their behavior. Thus, time dependent compactifications do not 
substantially modify the p-form spectrum; massless modes are still given by the de Rham 
cohomology of A4. 

C. The gravitational spectrum 

A key property of Kaluza-Klein reduction is that degrees of freedom in the full metric 
Gmn appear in lower dimensions as metric, vector, and scalar degrees of freedom. As 
in the p-form compactification discussed above, the masses of these fields depend on the 
properties of the compactification manifold A4. In contrast to the p-form case, the masses 
are not determined by the cohomology of A4, but by the existence of Killing fields on Ai. 
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The properties of Kaluza-Klein reduction, along with our discussion of chaos in Section 
ITT1 provide some useful simplifications. Since some of the metric degrees of freedom in the 
higher dimensional theory appear as zero- and one-forms, chaos arising from these degrees 
of freedom can be suppressed if they acquire a mass, just as in the conventional p-form case. 

As an explicit example of the reduction process, and of how chaos in higher and lower 
dimensions are related, we consider below the simple case with a single extra dimension jjjj]. 
The Kaluza-Klein reduction begins with a reparameterization of the metric, 

( e-^g^u + e 2 ^A^A u e 2 ^A„\ 

GMN ~ \ OA. o A. ' { '> 

where we assume that A^ and <fi are independent of the fifth dimension. We substitute 
this metric into the Einstein-Hilbert action and integrate over the fifth dimension. The 



coefficient q = a/2/3 is chosen so that the scalar field has a canonically normalized kinetic 
term in the resulting action, 

S = J R(g) - (<90) 2 - -^F^ sf^g d% (48) 

with F^ u = d IM A u — d u A fM . This describes Einstein gravity coupled to a scalar field (f>, and a 
vector field with ^-dependent coupling. It can be seen that all of the five-dimensional metric 
degrees of freedom in ()47|) are reproduced in this four-dimensional action. Furthermore, the 
vector term in the action possesses an exponential coupling to 0, of the type introduced 
in (PQ). Our starting point, the five dimensional pure gravity theory, is chaotic since the 
gravitational stability conditions cannot be satisfied for any choice of the Kasner exponents. 
After reduction, chaos also inevitably arises since the gravitational and one-form stability 
conditions cannot be satisfied simultaneously. Thus violations of the gravitational stability 
conditions in five dimensions can appear as violations of the p-form stability conditions in 
four dimensions. As we will discuss in more detail below, the preservation of chaos is not a 
generic feature of Kaluza-Klein reduction in dimensions greater than one. 

The example above is limited to a single extra dimension, and neglects metric modes 
that depend on the fifth coordinate. Below, we consider the general case, and calculate 
the effective masses of all Kaluza-Klein vector fields with an arbitrary number of extra 
dimensions. We will find that these masses are zero only when M. possesses Killing vectors. 



This calculation parallels standard treatments of Kaluza-Klein reduction 



231 ] . but in these 
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treatments the fact that Ai may not possess isometries is often not emphasized. For n > 1 
extra dimensions, we will generalize the decomposition (J4T)) using the vielbein formalism. 
We begin by defining one form fields e A = eu A dx M so that, 

ds 2 = e A e B VAB , (49) 

with t]ab the (4 + n)-dimensional Minkowski metric. The ej,/ are chosen so that, 

h^u = e^e/Vap, (50a) 

fmn &m C n 5 a b, (50b) 

e li a = K aj (x m )A? lt (x' i ), (50c) 
e m a = 0, (50d) 

where S a b the Euclidean flat space metric. The K a i are a basis for vector fields on Ai, indexed 
by j, that depend only on the compact coordinates x m . The coefficients in this expansion 
are the A? , which depend only the noncompact coordinates x^. The known as Kaluza- 
Klein vectors, will emerge after compactification as vector fields on the noncompact space 
R 1 ' 3 . The commutators of the K a i define a set of structure constants / J u 

[K j ,K k ] = f jk l K l . (51) 

The calculation is most conveniently carried out using an orthonormal basis {e a ,e a } given 
by, 

e a = e a e a = e a - K aj A j fl dx" , (52) 

in which the line element assumes the simple form ds 2 = e a e /3 r] ce j3 + e a e b S a b. In the event 
that some of the K a i are Killing fields on Ai, then the lower dimensional theory will possess 
a gauge symmetry. The Killing fields are generators of the isometry group of Ai, and 
this isometry group reemerges as the gauge group in the lower dimensional theory. This 
motivates the definition of a "field strength" F l as, 

F i = dA { + -p k A j A A k . (53) 
2 

In the general case, the Killing fields alone do not provide a full basis for vector fields on 
At. Thus, in addition to the massless modes (if any) of the gauge theory, there will also be 
an infinite set of massive gauge fields in the lower dimensional theory. 
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In order to derive the mass spectrum in the lower dimension explicitly, we may use the 
vielbeins to decompose the gravitational action in 4 + n dimensions. The spin connections 
are, 

L0 ab = LOab - \{V h K{ - V a K 3 h )Ay, (54a) 
= \KK^ a + V^,^, (54b) 

1 

2' 



= + 7^ K i F i^\ ( 54C ) 



where u a b and u a p are the spin connections defined by the metrics f mn on A4, and h^ v on 
R 1,3 , respectively Using these spin connections to compute the Ricci scalar, one obtains, 

R(G) = R{h) + R(f) - ^KiK^F^F^ - 2V \ c K j d) V {c K kd) A j "A*. (55) 

We see that a mass term for the A 7 has appeared. Upon integrating over the compact 
coordinates, one arrives at the Jordan frame action, 

S = J (w(f)R(h) - S(f) - \a jk F^F k ^ - feFAjj Vhd*x. (56) 

where, 

W(f) = [ yff^x (57a) 

JM 

S(f) = - [ R(f) y/ftTx (57b) 

JM 

a 3k = [ K ia K a k ^]d n x (57c) 

JM 

/3 jk = 2 [ V (c K ]d) V {c Kf ^fd n x. (57d) 

JM 

The W(f) factor may be removed by a rescaling of the metric h^ V) putting the action in the 
Einstein frame form. The S(f) term will yield a system of scalar fields, of which <p in our five 
dimensional reduction (|48j) is an example. Applying the Gram-Schmidt orthonormalization 
process to the K\ one can reduce a.j k = 5j k , giving the vectors a canonical kinetic term. 
Thus, we see Kaluza-Klein reduction of pure gravity results in a theory with scalars and 
vector fields, generalizing the n = 1 result discussed above. 

Of crucial importance to the present work is that massless Kaluza-Klein vectors are in 
one-to-one correspondence with Killing fields on A4 , or equivalently the zero eigenvalues of 
(3jk- This follows from the fact that j3j k functions as a mass matrix for the Kaluza-Klein 
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vector fields. Since (3jk is symmetric, we are guaranteed that m? will be real for all modes. 
In our discussion of p-form fields, we were able to apply powerful results regarding the 
Hodge-de Rham operator A that guaranteed that m 2 > 0, regardless of the topology and 
metric structure of M.. In the present situation, we have no guarantee that the masses of 
Kaluza-Klein vectors will satisfy m 2 > 0, or equivalently that the eigenvalues of the mass 
matrix are nonnegative. 

In this work we will assume that all eigenvalues of the mass matrix are nonnegative, so 
that m 2 > for all Kaluza-Klein vectors. In the general case, it is necessary to compute 
ctjk and /3jk for each manifold of interest, and then check that this assumption holds on a 
case-by-case basis. A simple example is provided by the n-torus T n . Realizing the torus 
as R n /Z n , with coordinates (6i . . . 9 n ) ranging on (0, 27r), a convenient basis for vector fields 
on T ra is provided by, 

, \ \/2f) a cosn#, for n > 0, 

K = ^ >< (58) 



(27t)"/ 2 



sinn#, for n < 



where n G Z, the 9 a are unit vectors associated to each coordinate, and (a,n) label each 
basis field, replacing the abstract indices used above. Substituting this into (JH7)) we find, 

®(a,n)(b,m) = S a b5 n m, (59a) 

(3(a,n)(b,m) = 2n 2 S ab S nm . (59b) 

The Kaluza-Klein vector fields are therefore canonically normalized, with masses \/2n for 
n = 0, 1, 2 ... , showing our assumption is valid in this case. More sophisticated examples 



may be found in the literature 23(. As we will explain in more detail below, the m 2 > 



assumption will enable us to treat p-forms and Kaluza-Klein vectors on the same footing. 



IV. SELECTION RULES FOR THE STABILITY CONDITIONS 

With the tools developed in the previous sections, we are now prepared to discuss condi- 
tions on Ai that result in controlled chaos. We will show that the gravitational, electric and 
magnetic stability conditions, introduced in Section m are modified by compactification. 
Not all of the stability conditions remain relevant, and only a subset need be satisfied to 
ensure that chaos is controlled. This subset is defined by the "selection rules" that are the 

24 



focus of this section. The selection rules that determine when a stability condition remains 
relevant are given for matter fields in Section liV Al and for gravitational modes in llVBI The 
selection rules, in turn, are determined by the de Rham cohomology (in the p-form case) 
and existence of Killing vectors (in gravitational case) of the compactification manifold M.. 
Here, we focus on discussing the origin of the selection rules. Once established, we will use 
them to find compactifications that control chaos in Section 

A. The p form selection rules 

The discussion in Section IIHI enables us to define selection rules for the electric and 
magnetic stability conditions. Each component of the p-form field results in an electric or 
magnetic stability condition, which expresses whether the energy density in that component 
scales rapidly enough to dominate the energy density of the universe and cause chaos. If this 
component gains a mass by compactification, then we have shown that it scales too slowly 
to be cosmologically relevant, and therefore we should ignore the corresponding stability 
condition. Thus we should ignore all electric and magnetic stability conditions involving 
indices that do not correspond to massless p-form modes. This results in the following 
selection rule, 

The p— form Selection Rule: When dim7^ s (M) = for some s, ignore the 
subset of p-form stability conditions, 

— T^-V<i> > (electric), (60a) 

v 

y^Pj - -Ajty < 1 (magnetic), (60b) 

with s Kasner exponents along the compact space A4. Retain only those stability 
conditions with s exponents along A4 and dim7i s (M.) ^ 0. 

In this section, we will use the results of previous sections to prove this rule, and give some 
simple examples of its use. 

This rule arises from considering the p-form modes that give rise to massless fields after 
compactification. A p-form A p gives rise to a massless r-form a r if and only if A p is of the 
form, 

A p = a r A[3 s (61) 
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where (3 S G TL S (A4), and r + s = p. Since d/3 s = 0, this gauge potential results in the field 
strength, 

F p+1 = (da) r+ i A P.. (62) 

When the energy density of this field is calculated, one finds a stability condition involving 
exactly s Kasner exponents along Ai, and the remainder along R 1,3 . Since only stability 
conditions of this type correspond to massless modes, they are the only ones that should be 
retained. 

The case (j62j) deals only with field strengths having at least one index along the noncom- 
pact space R 1,3 . In fact, the same selection rule applies when all indices of the field strength 
are along Ai. In this case, the field strength must satisfy the Bianchi identity and the Gauss 
law, 

dF p+1 = 0, SF P+1 = 0. (63) 

Field strengths of this type are commonly termed "nonzero modes" [^J . The conditions (jBHJ) 
imply that F p+ i is harmonic, and by Hodge's theorem the number of such forms is given by 
dim7^ p+1 (Ai). When dim7^ p+1 (Ai) vanishes, we cannot have a p-form field strength with 
all indices along Ai, and we should therefore delete the corresponding stability condition, 
with (p+ 1) Kasner exponents along Ai. Thus this case falls under the p-form selection rule 
as well. 

The selection rule may be illustrated by comparing compactification on a sphere S n and 
a torus T n . These manifolds encompass the best and worst case scenarios for controlling 
chaos through compactification. The sphere has the minimum number of massless modes for 
any orientable compact manifold, while the torus has massless modes for every dimension 
and involving every combination of indices on Ai. Therefore compactification on T n and 
S n will have very different influences on chaotic behavior. 

Compactification on T n does not modify any of the p-form stability conditions. The 
cohomology classes of the torus are, 

n ' 

dimft r (T n ) = „ " - , < r < n. (64) 
r\(n — r)\ 

If we realize the torus as R n /Z n , with coordinates 9\, . . . 9 n , then we may choose the following 
set of generators for the r th de Rham class, 

oj r = 69 h A d9 j2 A • ■ ■ A d6 jr , (65) 
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where {j r } are any set of distinct indices on T n . Therefore, massless modes exist for p-form 
fields with any combination of indices along the T n . Any p-form stability condition that 
appears in the noncompactified theory will remain in the compactified theory. 

By contrast, compactification on a sphere S n , with n > 1, deletes many of the stability 
conditions. The sphere has only two nonzero cohomology groups, each of unit dimension; 



The class TC° (S n ) is generated by the constant scalar function on S n , while the class H n (S n ) 
is generated by the volume form, 



where g mn is the metric, g = detg mn , and 8j the coordinates on the sphere. Massless modes 
therefore contain either no indices along the S n , or all indices at once. This implies that the 
only surviving stability conditions are those with either no internal Kasner exponents, or all 
internal Kasner exponents together. In the case where n > p, none of the internal Kasner 
exponents appear at all, and only those stability conditions involving Kasner exponents on 
R 1,3 survive. 

Our statement of the selection rule is the strongest possible in the generic case, and 
fortunately also the most conservative in terms of deleting the minimum number of stability 
conditions. Manifolds with a specific relationship between the frames a a appearing in (J2J) and 
the cohomology representatives of M may require that we delete additional p-form stability 
conditions. For example, consider the case in which M factors as M = Mi x M 2 , both 
topologically and metrically. A straightforward application of the selection rules results in 
retaining all stability conditions with s indices along M whenever dim7i s (M) 7^ 0. These 
stability conditions correspond to p-form modes with s x indices along Mi and s 2 indices 
along M 2 , with si + s 2 = s. However, in general a massless mode will not exist for every 
choice of si and s 2 , and therefore we may be able to delete additional stability conditions. We 
should only retain the even smaller subset of stability conditions with s± Kasner indices along 
Mi and s 2 indices along M 2 when dimH Sl (Mi) 7^ and dimH S2 (M 2 ) 7^ 0. Generally, 
however, we do not expect any special relationship between the aj and the cohomology 
classes. In this example, we have imposed the condition by hand that the aj point only 




for r = and r = n 



otherwise. 



(66) 



y/gd6i A d0 2 A • ■ ■ A d6, 



(67) 
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along exactly one of Ai\ or A\ 2 . Examples such as this one must be considered on a case- 
by-case basis, and lie beyond the scope of our selection rule. 

B. The gravitational selection rules 

The selection rules for the gravitational stability conditions arise in a manner similar 
to those for the p-form modes. We identify the degree of freedom corresponding to each 
stability condition, and then ignore the stability condition if the degree of freedom gains a 
mass through compactification. Unlike the p-form case, in general one must also check that 
the masses gained in this way satisfy m 2 > 0, as discussed at the end of Section Till CI In 
this way, we arrive at a selection rule for the gravitational stability conditions, 

The Gravitational Selection Rule: When M. possesses no Killing vectors, 
retain only the subset of gravitational stability conditions, 

Pi + Pj ~ Pk < 1) all triples i, j, k, (68) 

with all three Kasner exponents along R 1,3 , or all three along Ai, and ignore 
stability conditions with a mixture of exponents along both R 1 ' 3 and Ai. 

In proceeding, we are guided by the Kaluza-Klein reduced Jordan frame action (j^HJ). Clearly, 
gravitational stability conditions involving three Kasner exponents along the R ,3 should be 
retained, as the corresponding modes do not gain a mass from compactification. These 
represent the metric degrees of freedom in the lower dimensional theory. Stability condi- 
tions involving three Kasner exponents along the compact direction should also be retained. 
Physically, these correspond to metric degrees of freedom on compact space Ai. While these 
appear as scalar fields in the lower dimensional theory, they can result in a subtle form of 
chaos. Violations of these stability conditions appears as a chaotic system of interacting 
scalars in the lower dimensional theory. Thus, to ensure that all degrees of freedom are 
evolving smoothly to the big crunch, we should retain these stability conditions. 

Compactification can delete the "mixed" stability conditions, those with Kasner expo- 
nents along both R 1 ' 3 and Ai. These appear as the kinetic and mass terms for the Kaluza- 
Klein vectors in (|5fij) . When the compact space Ai possesses no Killing vectors, then these 
vector fields acquire a mass and become cosmologically irrelevant. When Ai possesses even 
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one Killing vector, then in general none of the mixed stability conditions can be discarded. 
This is because a single Killing field will generally involve all indices along Ai, and also 
results in a vector field with arbitrary indices on R 1 ' 3 . As in the p-form case, there can be 
special cases where additional stability conditions may be deleted. However, as we are more 
interested in the generic case we will not discuss examples here. 

V. EXAMPLES 

The previous sections have established that compactification allows us to ignore a number 
of the gravitational and p-form stability conditions. At this point, it is natural to ask if there 
are examples where enough stability conditions are deleted to control chaos. We show below 
that this is indeed the case, by giving explicit examples from both pure Einstein gravity 
and the low-energy bosonic sectors of string theory. We will discuss several solutions with 
controlled chaos in string models with M = 1 supersymmetry in ten dimensions, and will 
show how these solutions are interrelated by string duality relationships. 

For simplicity, we focus on regions on the Kasner sphere near what we will term "doubly 
isotropic" solutions. These are solutions in which a Kasner exponents take the value p a , and 
b take the value pb, with a + b = D — 1. In the absence of a dilaton, the Kasner conditions 
(J21) result in a quadratic equation for p a and pb, and therefore two solutions for each choice 
of a and b. When a dilaton is present, then there are two one-parameter families of p a 
and pb, which depend on the value of the dilaton "Kasner exponent" p^. Only models that 
are isotropic in three noncompact directions are of interest cosmologically, and so we fix 
a = 3. The Kaluza-Klein reduction of such models results in an isotropic universe with 
Pi . . .p 3 = 1/3, corresponding to a FRW universe dominated by a component with equation 
of state w = 1. 

It is important to emphasize that the specific examples that we will discuss are only 
representative points of an open region on the Kasner circle for which chaos is controlled, 
chosen so that the Kasner exponents assume a particularly simple, symmetric form. For a 
given compactification, the selection rules define a reduced set of stability conditions, which 
in turn define an open region of the Kasner circle for which all stability conditions are 
satisfied. When this open region is non-empty, then chaos is controlled. Thus, there will be 
choices of the Kasner exponents with controlled chaos in open neighborhoods of all of the 
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solutions discussed herein. 



A. Pure gravity models 

The simplest case, with n = 1 extra dimensions, is also a somewhat exceptional one. This 
is because there is exactly one compact one dimensional manifold, the circle S 1 . Regardless of 
the metric on the S 1 , it will always possess a Killing vector, and so no gravitational stability 
conditions can be deleted. Furthermore, Ti 1 (S 1 ) is nonzero, and so no p-form stability 
conditions are deleted. Therefore, all chaotic models remain chaotic when compactified on 
S 1 . To eliminate chaos when n — 1, we must consider a more general class of spaces than 
manifolds. 

For Einstein gravity without matter, a simple example that eliminates chaos when n = 1 
is given by the orbifold S 1 /Z 2 , previously discussed in Ref. Q]. If we take a coordinate 
9 on S 1 , ranging from [— 7r, it], then the orbifold results from identifying the S 1 under the 
reflection 9 — > —9. This takes G^g — > —G^g, and thus the Killing field is projected out, 
giving mass to all the Kaluza-Klein vectors. The resulting action for massless fields in four 
dimensions is then, 



in comparison with the classic Kaluza-Klein result (pfKJ) . Being only Einstein gravity with 
a scalar field, this theory is not chaotic. We will discuss compactification on this particular 
orbifold in more detail when we discuss string and M-theory solutions with controlled chaos. 

When we have n > 1 extra dimensions, then chaos can be eliminated by compactifying on 
a manifold without continuous isometries, and therefore without Killing vectors. This deletes 
the mixed gravitational stability conditions, as discussed in Section IIVBI The remaining 
stability conditions are always satisfied in the neighborhood of doubly isotropic solutions. 
This is subject to the assumption, discussed at the end of IIII C\ that the mass matrix for 
the Kaluza-Klein vector modes has no negative eigenvalues. 

While we have seen that the masses of Kaluza-Klein vectors are determined by isometric 
properties of Ai, there is a useful class of manifolds for which these properties are themselves 
determined by the topology, specifically by the de Rham cohomology. In this case, the 
gravitational and p-form selection rules are determined entirely by the cohomology of A4. 




(69) 
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These are the Einstein manifolds, for which, 

Rmn = ^9mni (70) 

with A arbitrary. When A4 is Einstein, the number of Killing vectors is given by dimli. 1 (A4) 
(3c]]. Chaos will thus be controlled in a neighborhood of doubly isotropic solutions when 
dim 7-^ {M) = 0. There are many examples of Einstein manifolds with this property; among 
them are the complex projective spaces CP" with the Fubini-Study metric, and the Calabi- 
Yau spaces. 



B. String models 

For string models with M = 1 supersymmetry in ten dimensions (Type I and heterotic), 
the simple class of doubly isotropic solutions is sufficient to give examples of solutions with 
controlled chaos. As we will discuss in more detail below, some of our solutions are related 
to others through standard string duality relationships. Interestingly, we find that theories 
with M = 2 supersymmetry (Type II) do not admit compactifications that lead to controlled 
chaos with doubly isotropic solutions. Unfortunately, we have found no examples where the 
compactification manifold could be a Calabi-Yau, although solutions with controlled chaos 
and Calabi-Yau compactification may exist for non-doubly isotropic choices of the Kasner 
exponents. 

In the following, we will always give the Kasner exponents in the Einstein conformal 
frame. Conventionally, the bosonic sector of string theory actions is presented in the "string 
frame" form, 

^t r in g = J (V 2 * [R(G^) + 4(«9$) 2 ] - £ e^Fl +1 j >/=GW d 10 x (71) 

where the A^ are the string frame couplings to the dilaton field, and the string frame 
metric. One arrives at the "Einstein frame" action by the transformation, 

ME) _ -*/2 G {S) (72) 

resulting in, 

Einstein = J ^(G**>) - {8<pf - ^ e^*ig + ij \Tg^ dPx (73) 
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where we have defined, 



(74a) 




8 - 2 Pj 



4 



) 



(74b) 



The field 0, which we shall refer to as the "dilaton" below, is canonically normalized, and 
the couplings between the dilaton and p-forms have transformed. In the following, we will 
always use the Einstein frame couplings, and so will drop the superscript (E) for clarity in 
notation. 

A theme common to our examples is the compatibility of our results concerning controlled 
chaos and the duality relationships connecting various string theories. In Section IV B II we 
first examine the Eg x Eg heterotic theory in detail. Through a combination of string duality 
relationships and compactifications we will be able to discuss its limits in eleven, ten, five 
and four dimensions. In particular, the S-duality relating the Eg x Eg heterotic string and 
M-theory 0, 0] is made apparent by relating the ten dimensional heterotic solution with 
controlled chaos and the compactification of eleven dimensional M-theory on S 1 /^- We 
will then discuss all compactifications of doubly isotropic solutions with controlled chaos for 
string theories with M = 1 supersymmetry in ten dimensions. We give four representative 
solutions, two each for the heterotic and Type I theories. We show that these four solutions 
organize into two pairs of solutions, related by the S-duality connecting the heterotic SO (32) 
and Type I strings 

It is important to keep in mind some features of the space of string solutions with con- 
trolled chaos. Each compactification we discuss, defined by the vanishing de Rham coho- 
mology classes, defines an open region on the Kasner circle where chaos is controlled. Our 
restriction to doubly isotropic models, in turn, takes a one dimensional "slice" out of this 
open region. In our examples, we give a representative point from the "slice" where the 
Kasner exponents assume a convenient and symmetric form. Thus we have found the com- 
pactifications that admit doubly isotropic solutions, but the choices of Kasner exponents are 
not unique. 
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theory 


spacetime 


dim 


Pi ■ --P3 


Pa ■ --P9 


Pio 


P<t> 


M-theory 


R 1 ' 3 xMx S7Z 2 


11 


-0.1206 


0.0662 


0.9644 




het E 8 x Eg 


R 1 ' 3 x M 


10 





1/6 






"braneworld" 


R 1 3 x SVZ 2 


5 


0.0105 




0.9686 


0.2486 


FRW 


R 1 ' 3 


4 


1/3 






a/273 



TABLE I: The solution discussed in the text for the Eg x Eg heterotic string in various dimensions. 
The Kasner exponent is a combination of the dilaton and volume modulus for the compactified 
dimensions. Chaos is controlled provided that H? (Ai) = 0. 

1. The heterotic string and M-theory 

Here we focus on the Eg x Eg heterotic theory, in the neighborhood of a specific choice of 
Kasner exponents. Using string duality relationships and compactification, we will discuss 
the various guises of this solution in eleven, ten, five and four dimensions, summarized in 
Table HI While the solution we discuss also controls chaos for the SO (32) heterotic theory in 
ten dimensions, string dualities for this theory do not enable us to discuss the five dimensional 
and M-theory limits. To begin, we consider the heterotic theory in Einstein frame, where 
it contains the metric Gmn, dilaton 0, one-form A\ and two-form B 2 . The dilaton couples 
to the one- and two-forms via exponential couplings of the type with Ai = — l/v2, 
and A 2 = Before compactification, one finds violations of the electric and magnetic 

stability conditions for the choice of Kasner exponents, 

Pl ...p 3 = 0, p4...P9 = l/6, Pt = V$/Q (10D). (75) 

We assume that p\ . . .p 3 lie along the noncompact spacetime R 1,3 , and p 4 . . .p 9 lie along the 
compact manifold Ai. In this solution, the magnetic stability conditions are violated for the 
magnetic component of H3 = di?2 with all three indices along Ai. None of the gravitational 
stability conditions are violated. 

Applying the selection rules introduced in Section IIVI we find that chaos is controlled 
by compactifying on a six-manifold Ai with H? (Ai) = 0, such as S 6 or CP 3 . The choice 
of CP 3 has the advantage that it has 7i l (CP 3 ) = 0, and is an Einstein manifold if given 
the Fubini-Study metric. This manifold may therefore be a useful starting point for models 
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that differ from the doubly isotropic ones. Unfortunately, a Calabi-Yau space will always 
have dim7i 3 (A4) 7^ 0, and is thus unsuitable for rendering this solution non chaotic. 

The four dimensional limit of the solution (|75|) possesses a simple form. This is obtained 
by compactifying on the six-manifold Ai, resulting in, 

Pi ...P3 = 1/3, (4D). (76) 

This describes a collapsing, flat FRW universe dominated by a perfect fluid with w = 1. 
The w — 1 component is a combination of the dilaton and the volume modulus arising from 
the Kaluza-Klein reduction of the heterotic theory on the six-manifold A4. 



String duality relationships imply that the (strongly coupled) 



x Eg heterotic theory 



is obtained by compactifying M-theory on the orbifold S 1 /^ 3 it l32t l33f| . Phenomenology 
implies that the orbifold S 1 /Z 2 is somewhat larger than the compactification six manifold 
M.. Thus, depending on the scale of interest, the strongly coupled heterotic theory can 
appear four, five, or eleven dimensional. The five and eleven dimensional limits of the 
solution (|75j) are not as simple as the ten and four dimensional views, but are nonetheless 
instructive. 

This duality relationship implies that the heterotic string in ten dimensions can be de- 
scribed by eleven dimensional M-theory on R 1 ' 3 x M. x S 1 /^. The eleven dimensional 
lifting of (J7H|) to M-theory yields Kasner exponents whose precise expression is not very 
illuminating, but whose approximate numerical values are, 

Pi...P3 = -0.1206, P4---P9 = 0.0662, p 10 = 0.96442 (11D). (77) 

This describes a rapidly shrinking orbifold, a slowly contracting six-manifold Ai, and a 
slowly expanding noncompact space. 

To see that the compactification of M-theory on R 1,3 x Ai x S 1 /Z 2 leads to controlled 
chaos requires us to consider some subtle features of the theory. The bosonic sector of M- 
theory includes only the graviton and a four form field strength F4 = (IA3. Our selection 
rules are only strictly applicable to the case where the compactification space is a manifold, 



an assumption which fails to include or 



here follows usual techniques 0, for determining the spectrum after compactification. 

Specifically, we compactify M-theory on S 1 , and then impose the identification 9 — > —9. 
The presence of the Chern-Simons term A 3 A F 4 A F 4 in the M- theory Lagrangian requires 
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Difolds such as S 1 /Z 2 . The approach most convenient 



that F4 — > — F4 under parity transformations, of which the identification 6 — > —0 is an 
example. Thus, the massless components of F4 on R 1 ' 3 x M. x S 1 /^ are those with exactly 
one index along the S 1 /Z 2 . 

Before compactification, the eleven-dimensional solution given above violates the p-form 
stability conditions for three components of the four form field. The first and second are 
electric, with the first having three indices along R 1 ' 3 , and the second having two along R 1,3 
and one along A4. The third is magnetic, with one index along S 1 /^ and three along M.. 
The magnetic component is rendered massive by the condition H? (M) = 0, and so we can 
neglect this stability condition. The two electric components are rendered massive since 
they do not have exactly one index along the S 1 /Z 2 , and their stability conditions can be 
neglected as well. Thus, chaos is controlled in the M-theory limit. 

The five dimensional guise of our solution, obtained by Kaluza-Klein reducing the eleven 
dimensional form on the six- manifold A4, describes a "braneworld" with structure R 1 ' 3 x 
S 1 /Z 2 . This yields the solution, 

Pl ...p 3 = 0.01048, p 10 = 0.9686, p^ = 0.24804. (5D) (78) 

The scalar field ip is the volume modulus of the six-manifold M.. This solution describes 
a nearly static R 1,3 and a rapidly contracting orbifold. This solution bears a suggestive 
similarity to the set-up studied in the ekpyrotic/cyclic scenario. In these models, near the 
big crunch, the five-dimensional spacetime approaches the Milne solution, 

Pi...p 9 = 0, pio = 1. (79) 

This solution is in fact on the boundary of the open region of the Kasner sphere for which our 
example solution (|75j) exhibits controlled chaos. This is predicated on the assumptions that 
w = 1 in the four-dimensional theory all the way to the big crunch, and that compactification 
is the only mechanism for controlling chaos. In the ekpyrotic/cyclic scenarios, there is a 
long w ^> 1 phase during the contraction, in which the energy density in p-form modes 
is exponentially suppressed l^j]. This suppression further reduces the time td om at which 
dangerous modes can formally dominate the universe. Therefore, in the full model, the onset 
of chaos will be delayed far beyond what our estimates, based only on the compactification 
mechanism, would suggest. 
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sol'n 


Pi ■ --P3 


PA---P9 


P<t> 


theories 


zero betti 


A 





1/6 


a/5/6 


heterotic 




B 





1/6 


-y/KJE 


Type I 


h 


C 


1/3 





a/2/3 


Type I 


b\M 


D 


1/3 





-a/273 


heterotic 


b\M 



TABLE II: Representative string theory solutions with controlled chaos and isotropic behavior 
along the noncompact and compact directions. Each compactification leads to open regions of the 
Kasner circle with controlled chaos, and we have given a representative point for each open region 
here. The string theories which exhibit controlled chaos for each solution are shown, as well as the 
Betti numbers bj = dimTC J {M) of M that are required to vanish. 

2. The heterotic and Type I strings 

Having focused in detail on a single compactification of a single string theory, we now focus 
on finding all compactifications with controlled chaos and doubly isotropic Kasner exponents. 
There are four doubly isotropic examples with controlled chaos, with representative choices 
of the Kasner exponents summarized in Table ILT1 The E 8 x E 8 and SO (32) heterotic theories 
exhibit the same chaotic behavior, since their p-form spectrum and couplings to the dilaton 
are identical; these theories differ only in the gauge groups for their non-abelian gauge 
multiplets. One may also include the ten-dimensional (noncritical) bosonic string, which 
contains only the Neveu-Schwarz fields of the heterotic string and no gauge fields. One finds 
that chaos is controlled in the ten dimensional bosonic string in the same solutions (A and 
D) as in the heterotic string. 

In the absence of any compactification, these examples are all chaotic. None of them 
suffer from gravitational chaos, and in all cases the chaotic behavior arises from the p-form 
fields alone. Upon compactification to four dimensions, these models all result in a FRW 
universe dominated by a free scalar field with w = 1. 

The examples given in Table [H] include not only models that go to weak coupling at the 
crunch, (A and C) but also models where the dilaton runs to strong coupling (B and D). 
The fact that the solutions include both those where the string theory goes to strong and 
weak coupling is interesting from a model building perspective. The dilaton is never static 
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in the solutions discussed here, a feature also found in some cosmological models based on 
string theory. In the ekpyrotic/cyclic models, for example, the string coupling goes to zero 
at the big crunch. In pre big-bang models, on the other hand, the dilaton goes to strong 
coupling at the crunch. Thus the controlled chaos mechanism may be relevant to both 
scenarios. 

The heterotic SO (32) and Type I theories are related by an S-duality transformation, 
and this symmetry is respected by our examples here. Under this duality, the string frame 



actions of the heterotic SO (32) and Type I theories are related by |25(, 

Gmn = e-**G<**> (80a) 

$C0 = (80b) 

with the p-form fields remaining unchanged. Carefully working through the resulting trans- 
formation of the Einstein frame Kasner exponents, one finds that the spatial Kasner ex- 
ponents are unchanged, while p^ = —p^^- Therefore, ^-duality exchanges the pairs of 
solutions A <-> B and C <-> D. 

The properties of string theories regarding controlled chaos appear correlated to their 
supersymmetry properties in ten dimensions. The Af = 1 theories, (heterotic, Type I, and 
M-theory on S 1 /^) possess simple compactifications that control chaos. The Type IIA/B 
theories and uncompactified M-theory, with Af = 2 supersymmetry, have no doubly isotropic 
solutions with controlled chaos. As we have not exhaustively examined the Kasner sphere, 
we cannot say for certain whether there exist solutions that control chaos for the Af = 2 
theories. 

It is natural to expect that the Af = 1 and Af = 2 string theories will have different 
characteristics with respect to controlled chaos. There is a useful formulation of the dynamics 
of gravity near a big crunch, discussed briefly in Section|nJ In this formulation, the dynamics 
of metric and p- form fields is recast as the motion of a billiard ball in a hyperbolic space, 
undergoing reflections from a set of walls. The walls correspond to p-form kinetic terms and 
curvature terms in the Einstein equations. The positions and orientations of these walls are 
identical for all of the Af = 1 theories, and different from the common set of walls shared 
by the Af = 2 theories |6j. Our suppression of the energy density in massive p-form and 
gravitational modes amounts to "pushing back" these walls. Thus, it is not surprising that 
we should find that Af = 1 and Af = 2 models have different characteristics with respect to 
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controlling chaos. 



VI. CONCLUSIONS 

The results presented here build on the many years of previous research in the behavior 
of general relativity near a big crunch. Previous research has primarily focused on "local" 
properties of theories with gravity, such as the dimensionality of spacetime, or the types and 
interactions of matter fields, and has revealed how these influence the emergence of chaos. 
Here we have investigated "global" features, in particular the topology of spacetime. We 
have found that these features can lead to a suppression of chaos in many models of interest. 
The control of chaos can be expressed simply in terms of selection rules for the gravitational 
and p-form stability conditions. These in turn can be used to find compactifications of 
chaotic theories in which chaos is suppressed right up to the quantum gravity regime. 

Our results bear an intriguing connection to some cosmological models that are founded 
on current ideas in string and M-theory. Among the simple examples of string theory 
solutions with controlled chaos, we find those that resemble both the ekpyrotic/cyclic and 
pre-big bang scenarios. For future models, this work suggests a method to control chaotic 
behavior near a big crunch that does not require postulating additional interactions and 
matter fields, or depending on higher order corrections to the Einstein equations. While 
this work sheds no light on the behavior of these models in the quantum gravity regime or 
through the big crunch/big bang transition, it provides a natural mechanism that ensures 
that the universe evolves smoothly so long as classical physics may be trusted. 

Recent work suggests that maintaining this smooth contraction during the classical regime 
may be sufficient to allow a nonsingular quantum evolution through a big crunch/big bang 
transition. One approach to this problem begins from the fact that, in string and M- 
theory, the degrees of freedom during the quantum regime are very different from those of 
the classical regime studied here. The fundamental degrees of freedom are extended objects, 
such as strings and branes. As one approaches the scale set by their tension, classical general 
relativity breaks down, and these extended objects become the relevant degrees of freedom. 
In particular, it is the evolution of these strings and branes that one should study near the 
big crunch. Working within the context of the ekpyrotic/cyclic scenario, it was found in 
Ref [2Cj that if the universe is sufficiently smooth and homogeneous at the beginning of 
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the quantum regime, the fundamental excitations (M2 branes) evolve smoothly through the 
big crunch with negligible backreaction. This suggests that a sufficiently smooth "in" state 
can evolve through the big crunch to a smooth "out" state, precisely what one requires for 
cosmology. This result complements the present work. The mechanism described herein 
can be viewed as providing the required conditions for smooth classical evolution before the 
Einstein equations break down, preparing the universe for nonsingular quantum evolution 
through the big crunch. 

Our results have further implications for high energy theory and phenomenology. String 
models and M-theory require compactification in order to produce the correct number of 
observed noncompact dimensions. Obtaining the correct low energy physics, such as M = 
1 supersymmetry in four dimensions or the correct number of lepton generations, puts 
constraints on the compactification manifold A4, many of which are topological in nature. 
Controlling chaos through compactification in cosmological models with a collapsing phase 
places additional constraints on M. We are currently investigating whether these two set 
of constraints are compatible. For example, the existence of solutions with compactification 
on a Calabi-Yau space would suggest that chaos can be controlled in string models with a 
realistic low energy spectrum. 

These results also inspire more speculative scenarios. When the universe enters a chaotic 
regime, the Kasner exponents will undergo an infinite number of "jumps" to different points 
on the Kasner sphere as the big crunch is approached. We also might expect that the topol- 
ogy of M, is changing at the same time. For example, there are situations in string theory 
where the topology of M. can change dynamically, such as the conifold or flop transitions. 
If the combination of Kasner exponents and topology lead to controlled chaos, then the 
universe will subsequently contract smoothly to the big crunch. In this way, the universe 
will have dynamically selected not only some properties of M., but also a "preferred" cos- 
mological solution near the big crunch. Analysis of such a scenario would require a much 
deeper understanding of cosmology in the quantum gravity regime than is currently avail- 
able, clearly an important topic for further research. 
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Appendix 1: Hamiltonian Formulation of p form Dynamics 

Below, we treat the case of a general p-form field with coupling to the dilaton <f), in 
a collapsing universe. We concern ourselves with the case where spacetime is isotropic, 
as this is the cosmologically relevant situation after compactification. We show that far 
from the big crunch, the energy density in massive fields evolves like that of a pressureless 
fluid, p ~ 1 / (comoving volume) . We will first recast the p-form dynamics in Hamiltonian 
form. This allows us to apply the virial theorem and stress energy conservation to obtain 
the scaling in energy density far from the crunch. The p-form action with mass term and 
dilaton coupling is, 



S 



- , 1 / (dA-dA + m 2 A-A)e X(t, V^Gd D x, (81) 

where we have fixed the coordinate gauge so that ds 2 = —n 2 dt 2 + r )^dx^dx k . We choose 
the canonical coordinates to be the gauge potential A ai ^ ap . The corresponding canonical 
momenta are, 

U ji...j p = _ F tji-j Pne \<f>^ ( 82 ) 
Passing to the Hamiltonian, we find, 

H= l -Jfi{n-U + 7 e 2A *F^ • + ie 2 ^m 2 A ■ A - ie 2 ^A ta2 ... ap d a U aa ^) d D ~ x x 

(83) 

where we use a rescaled lapse function n = ne~ x< ^/ */7, and denote the magnetic components 
of F p by F^ B \ Dot products are taken with respect to the metric Jmn- The last term in the 
integral shows that the "electric" gauge field modes A t j 2 _j v _ x appear as Lagrange multipliers 



necessary to enforce the Gauss's law constraint, but are otherwise nondynamical [34J. We 
will choose the Coloumb gauge, in which A t j 2 ^j p = 0, and drop this constraint term from 
now on. 
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The Hamiltonian is in fact exactly that of a set of simple harmonic oscillators. After 
decomposing the functions A p in an appropriate orthonormal set of Fourier components, 
different Fourier modes decouple and the Hamiltonian is quadratic in II and A. The electric 
field modes appear as the kinetic terms, and the magnetic field and mass terms correspond 
to the potential of the oscillators. The oscillator potential is time dependent, both due to 
the appearance of e 2A<?i 7 and individual metric components in the magnetic and mass terms. 

We are primarily interested in the dynamics of the p-form far from the big crunch. 
In this regime, we may view the changing scale factors as slowly varying parameters in 
our Hamiltonian. They will change the spring constants on a timescale given by t, the 
proper time to the big crunch. The dynamical timescale (typical period) for the oscillator 
Hamiltonian is given by the mass term. Thus, we expect that the fractional change in the 
fundamental frequencies to of the oscillator system over a typical cycle will be, 

5uj 1 . , 

84 

uj ujt K ' 

Provided we are at a time t ^> we will be in the adiabatic regime, and we may take 
the oscillation frequencies to be constants. This corresponds to the m/H 3> 1 regime that 
we have concerned ourselves with throughout this paper. 

To find how the energy density scales with time, we apply the virial theorem and stress 
energy conservation. The adiabatic condition ()84j) implies that we may neglect the time 
variation of the metric and dilaton over a single cycle. For our Hamiltonian, the virial 
theorem then implies that the time average (•) of the potential energy is equal to that of 
the kinetic energy, or 

(n • n> = e 2X ^(F^ ■ F< s > + m 2 A ■ A). (85) 

In the virialized system, there are two possible regimes, corresponding to either the -F^ 
term or the m 2 A ■ A term dominating. We will consider both of these cases in turn. 
The stress energy for the p-form field is, 



p \4> 

rp V _ ° 



(P+ l)F, a2 ... ap F^-^ - +pm 2 A lia2 ... ap „ 1 A^- a ^ - —8J 'A 2 



(86) 



(p+1) 

We will find it convenient to break the stress energy into three parts, 

r; = i< B ); + 2<V + T (B) ; (87) 
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corresponding to the energy in electric modes, magnetic modes, and the mass term. It is 
sufficient to consider a single component of F p , since different components will be uncor- 
related and therefore will have vanishing time average. The electric modes give rise to a 
contribution, 



3 A0 



t 1 — o/„ i 1 M 1-^(^)1 X 



2(p + l)!' 



while the magnetic modes give, 



D A</> 



2(p + l)! 



F (B) >< 



and the mass term yields, 



M 2(p+l)! 



m 2 A 2 x 



■I if has index /i, 

-1 otherwise. 

+1 if has index /i, 

- 1 otherwise. 

-1 if A has index p,, 

— I otherwise. 



(88) 



(89) 



(90) 



Note that F^ cannot have any timelike indices, nor can A thanks to our gauge choice. 
Thus the contributions to the energy density p = — To° are all positive. 

First, we consider the case where the mass term dominates in the virial relationship 
()85|) . This corresponds to inhomogeneities in the p-form field being negligible. The virial 
result implies that (\F^\) = m 2 (A 2 ). Due to our gauge choice, A and F^ E > have the same 
combination of p spatial indices, and therefore contributions to the pressure components 
7y coming from T^/ and T^/ exactly cancel. The vanishing pressure reveals that the 
effective equation of state is that of dust, w = 0. 

When the magentic terms dominate the virial result 1)85)1 . we obtain a slightly different ef- 
fective equation of state. Here it is necessary to average over polarizations of F p , since unlike 
F^ and A, F^ and F^ do not enjoy any relationships between their indices. Regardless 
of polarization, the sum of stress energy tensors and has vanishing trace. This, 
combined with isotropy, implies the pressure components are given by = —Tq°/(D — 1). 
This corresponds to the equation of state of radiation, which in four dimensional spacetime 
is w = 1/3. 

Physically, this result may be understood in simple terms. Far from the big crunch, the 
contraction of space is very slow in comparison to the mass of the p-form field. Thus, the 
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corresponding particles are far from being relativistic, and behave as a dust of approximately 
comoving mass points. Their energy density therefore scales in inverse proportion to the 
comoving volume. The case where magnetic components dominate the virial relationship 
corresponds to a relativistic gas of particles. This yields the equation of state of radiation, 
as we expect. 
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